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^ \ ABSTRACT 

O [ Aims. We investigate the role of convection in the formation of atomic and molecular lines in the atmosphere of a red giant star. For 

this purpose we study the formation properties of spectral lines that belong to a number of astrophysically important tracer elements, 
including neutral and singly ionized atoms (Lii, Ni, Oi, Nai, Mgi, All, Sii, Sin, Si, Ki, Cai, Can, Til, Tin, Cri, Cm, Mm, Fei, 
Fen, Co i, Ni i, Zni, Srn, Ban, and Eun), and molecules (CH, CO, C 2 , NH, CN, and OH). 

Methods. We focus our investigation on a prototypical red giant located close to the red giant branch (RGB) tip (r e fj=3660 K, 
logg=1.0, [M/H]=0.0). We used two types of model atmospheres, 3D hydrodynamical and classical ID, calculated with the C0 5 B0LD 
and LHD stellar atmosphere codes, respectively. Both codes share the same atmospheric parameters, chemical composition, equation 
of state, and opacities, which allowed us to make a strictly differential comparison between the line formation properties predicted in 
3D and ID. The influence of convection on the spectral line formation was assessed with the aid of 3D-1D abundance corrections, 
£f>^ ' which measure the difference between the abundances of chemical species derived with the 3D hydrodynamical and ID classical 

model atmospheres. 

Results. We find that convection plays a significant role in the spectral line formation in this particular red giant. The derived 3D-1D 
^ — ■ abundance corrections rarely exceed ±0.1 dex when lines of neutral atoms and molecules are considered, which is in line with the 

previous findings for solar-metallicity red giants located on the lower RGB. The situation is different with lines that belong to ionized 
atoms, or to neutral atoms with high ionization potential. In both cases, the corrections for high-excitation lines (x > 8 eV) may 

amount to A 1D _i D 0.4 dex. The 3D-1D abundance corrections generally show a significant wavelength dependence; in most cases 

they are smaller in the near-infrared, at 1600-2500 nm. 
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S— i | 1 . Introduction to which extent convection and other nonstationary phenomena 

may influence the observable properties of red giant stars. 
Current three-dimensional (3D) hydrodynamical codes have v tQ ^ Qnl a few smdies haye focused Qn ^ inyes _ 

brought a new level of realism to the modeling of stellar at- .... e u C. e a ■ » - A m u a 

& , ... ., , • rf r tigation of observable properties of red giants with 3D hydro- 

mospheres, making it possible to assess the influence of van- d ical stellar model atmospheres . In one of the early at- 
ous nonstationary phenomena on the observable properties of Kucinskas et al. (2005) have found that convection may 

various classes of stars. Since red giants are amongst the bright- .. , , u tU , t e tU A ■ * *™ 

.. ,. ,. • noticeably alter the thermal structure of the red giant atmo- 

est objects in intermediate age and old populations, precise un- here whichi ^ affect ^ ^ 

derstanding of their internal structures and observable proper- tion of ^ eme in radiation field xhis ield a differ ence 

ties is of key importance tor studies of stellar populations in the rnonQ -ujuju* / i a- » a 

_ , ' r . _ ,. . , , , , • , „ of ~ 0.2 - 0.3 mag in broad-band photometric colors predicted 

Galaxy and beyond. Three-dimensiona hydrodynamical stellar fe ^ 3D h drodynamical and classical 1D models . It was 

model atmospheres may prove especially useful in this context ^ shown fe CoUet gt aJ ( 2QQ1) ^ convection in the atmo . 

as they can provide important theoretical insights about how and , c „ , . , • . /T , Annt\ cmni^ 

j f f & spheres of somewhat warmer red giants (T e g « 4700 - 5100K, 

logg - 2.2, [Fe/H] = 0.0 to -3.0) may lead to significant 
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ities. Remarkably, at [M/H] = -3.0 this may result in differences 
of up to -1.0 dex in the abundances of chemical elements de- 
rived with the 3D hydrodynamical and classi cal ID model atmo- 
sphere s. T his result was late r conf irmed by iDobrovolskas et all 
(120101) and llvanauskas et all (1201 Ol) , who used different 3D hy- 
drodynamical stellar atmosphere and spectrum synthesis codes, 
C0 5 B0LD and Linfor3D, respectively, together with the atmo- 
spheric parameter s of red giants similar to those utilized by 
IColiet et al.l (|2007). All these studies unequivocally point to the 
fact that the role of convection in the atmospheres of red gi- 
ant stars may be significantly more important than previously 
thought, especially at metallicities [M/H] < -1.0. 

Aiming at a more systematic study, we have recently started 
a project to investigate internal structures and spectral properties 
of red giants across the HR diagram with the aid of 3D hydro- 
dynamical model atmospheres. For this purpose we used the 3D 
hydrodynamical stellar atmosphere code CO 5 BOLD to produce a 
grid of red giant model atmospheres, which are now available 
as part of the CIFIS T grid of C0 5 B0LD 3D model atmospheres 
(ILudwig etalJ [2009. the grid is constant ly updated with new 
models). In the first study of this series (ILudwig & Kucinskas! 
1201 2|) we investigated physical properties of convection in the 
atmosphere of a red giant located close to the red giant branch 
(RGB) tip (T eff ~ 3660 K, logg = 1.0, [M/H] = 0). In the 
present study we extend our work on the same red giant and fo- 
cus on the influence of convection on the formation of various 
atomic and molecular lines. Our main goal is to make a detailed 
comparison between the line strengths predicted by the 3D hy- 
drodynamical and classical ID models, and discuss the conse- 
quences for stellar abundance studies. General properties of the 
spectral energy distribution and photometric colors of this partic- 
ular red giant will be discussed in a companion paper (Kucinskas 
et al. 2012, in preparation). 

The paper is organized as follows. The model atmospheres 
used in this work are described in Sect.[2j where we also outline 
the details of spectrum synthesis calculations and the concept 
of 3D abundance corrections. The main results are presented in 
Sect. [3] including a brief discussion of the basic properties of 
spectral line formation in the presence of convection, and the dif- 
ferences between the predictions of the 3D hydrodynamical and 
classical ID models. The conclusions are presented in Sect. [4] A 
more detailed analysis of the abundance corrections derived in 
this work as well as additional background information is given 
in Appendices A, B, and C. 

2. Stellar model atmospheres and spectral line 
synthesis 

2.1. Model atmospheres 

We used red giant model atmospheres calculated with the 3D 
hydrodynamical C0 5 B0LD and ID stationary LHD model atmo- 
sphere codes. Both models were computed for the same atmo- 
spheric parameters (r eff =3660K, log g= 1.0, [M/H]=0.0) and 
elemental composition, and with the same equation of state 
and opacities; they also share the same numerical radiative 
transfer scheme. Additionally, we also used an average (3D) 
model, which was produced by spatially averaging the 3D 
model structures at constant Rosseland optical depth. The (3D) 
model does not contain explicit information about the horizon- 
tal fluctuations of thermodynamical and hydrodynamical quan- 
tities (e.g. temperature, pressure, velocity) that are present in 
the full 3D model atmosphere. Therefore, the comparison of 
the predictions of the full 3D and average (3D) models allows 



one to estimate the relative importance of the horizontal tem- 
perature fluctuatio ns on the spectral line formation (see, e.g. 
ICaffau et al.l 1201 ll for a more detailed discussion). Since the 
thr ee types of models used in this work are described in detail 
by ILudwig & Kucinskas! d2012l) . we only briefly summarize the 
most essential aspects of their calculation in the following sub- 
sections. 

2.1.1. Three-dimensional hydrodynamical model and 3D 
snapshot selection 

The red giant model atmosphere used in this work was cal- 
culated with the 3D radiation hydrodynamics code C0 5 B0LD 
dFrevtag et al.ll2012l) . The C0 5 B0LD code uses a Riemann solver 
of Roe type to calculate the time evolution of the hydrody- 
namical flow and the radiation field on a 3D Cartesian grid. 
The model was computed using a grid of 150 x 150 x 151 
mesh points (ixyx z), which corresponds to a physical box 
size of 15.6 x 15.6 x 8.6 Gm 3 . The radiative transfer is based 
on monoc hromatic opacities from the MARCS stellar atmosphere 
package (Gustafsson et al. 2008) which, to speed up the cal- 
culations, were grouped into five opa city bins (for details on 
the opacity binning approach see, e.g. i Nordlundl l 19821 : ILudwig 
119921: ILudwig et alJfl99lt IV6glerll2004l) . Solar elemental abun- 
dances as given by iGrevesse & Sauval (Il998h were assumed, 
except for carbon, nitrogen, and oxygen, for which the follow- 
ing values were us ed: A(C)=8.41, A(N)=7.8, and A(0)=8.67 
(Caffauet al. 2008). The model calculations were made assum- 
ing local thermodynamic equilibrium, LTE (for mo re details 
about the model setup see Ludwig & Kucinskas 2012). 

After the initial relaxation to a quasi-stationary state, we ran 
the model simulations to cover a span of ~6 x 10 6 sec (~70 days) 
in stellar time. This corresponds to approximately seven convec- 
tive turnover time£] in the atmosphere of this particular red giant. 
From this sequence of relaxed models we selected 14 3D model 
structures computed at different instances in time (snapshots). 
Individual snapshots of this 14-snapshot subset were then used 
to produce average (3D) model and spectrum synthesis calcu- 
lations (Sect. 12. lT2l and l2~2l respectively). Snapshots of this 14- 
snapshot ensemble are separated by 5 x 10 5 sec (~6 days) in stel- 
lar time, which allows one to assume that in this subset they are 
statistically uncorrected. The snapshots were selected in such 
a way that most important statistical properties of the snapshot 
ensemble, such as the average effective temperature and its stan- 
dard deviation, mean velocity at optical depth unity, mean veloc- 
ity profile and residual mass flux profile, would match those of 
the entire 3D model run as closely as possible. 

2.1.2. Average (3D) model 

For each of the selected snapshots, the corresponding average 
3D model, (3D), was obtained by spatially averaging the thermal 
structure of the 3D model box on surfaces of equal Rosseland 
optical depth. We averaged the fourth moment of temperature 
and first moment of gas pressure to preserve the radiative prop- 
erties of the original 3D mo del, according to the prescription 
given in Steffen et al. ( 1995). The product of this procedure is 
a ID model atmosphere that retains the averaged vertical pro- 
file of the thermodynamical structure of the 3D model, but lacks 
explicit information about the horizontal inhomogeneities. This 



1 as measured by the Brunt- Vaisala and/or advection timescales, the 
latter equal to the time needed for the convective mat erial to cross 1.8 
pressure scale heights (see Ludwie & Kucinskas 2012). 
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was done for the 14 3D model snapshots selected at different in- 
stants in time (Sect. l2.lTTI ), thus obtaining a sequence of 14 (3D) 
models. 



2.1.3. One-dimensional (1D) model 

The ID hydrostatical model was calculated with the LHD code, 
using the same atmospheric parameters, elemental abundances, 
opacities, and equation of s tate as in the 3D model calculations 
described above (see, e.g.. ICaffau et al.l [20081 for more details 
on the LHD cod e). Convection in the LHD models was treated ac- 
cording to the Mihalas (fl978h formulation of the mixing-length 
theory. The LHD models were calculated for the mixing-length 
parameters q-mlt = 1.0 and 2.0, to investigate the influence of 
Q"mlt on the properties of line formation in the ID models (see 
Sect.P 



2.2. Three- and one-dimensional spectrum synthesis 
calculations 

2.2.1. Chemical elements and spectral line parameters 

We used fictitious spectral lines, i.e., lines of a particular chem- 
ical element (molecule) for which the central wavelength, A c , 
the excitation potential of the lower level, and line equiva- 
lent width, W, were selected arbitrarily. This approach allowed 
us to cover a range in A c and x to quantify the trends of line 
formation properties in the 3D hydrodynamical and ID classi- 
cal model atmospheres with respect to these two line parame- 
te rs. Conceptually, this metho d can be traced back to the work 
of ISteffen & Holwegeri 12 002) and was applied in several later 
studies too (e.g ICollet et al] 12007. iDobrovolskas et al.l 12010. 
llvanauskas et al.ll2010t) . 

Synthetic line profiles were calculated for the following 
tracer species of astrophysical interest: 



- neutral atoms: Li i, N i, O i, Na i, Mg i, Al i, Si i, S i, K i, Ca i, 
Ti i, Cr i, Mn i, Fe i, Co i, Ni i, and Zn i; 

- ionized atoms: Sin, Can, Tin, Cm, Fen, Srn, Ban, and 
Eun; 

- molecules: CH, CO, C 2 , NH, CN, and OH. 

Fictitious lines were calculated at three wavelengths, A c = 
400, 850, and 1600 nm. The first two were selected to 
bracket the wavelength range accessible with modern high- 
resolution optical spectrographs (e.g., UVES/GIRAFFE@VLT, 
HARPS @ES03.6m, HIRES @Keck), whereas the third cor- 
responds to the wavelength reachable with similar instru- 
ments in the near-infrared //-band (e.g., CRIRES@VLT, 
NIRSPEC@Keck). In the case of molecules, however, we used 
the wavelengths of real molecular bands in the blue part of the 
spect rum, whenever avai lable: the A-X OH transition at 315 nm 
(e.g.. iBessell et alJl2004 . the A-X transition of NH at 336 nm 
(e.g.. ISpite et all I2005TT the B-X transition of CN at 388 nm 
(lLambert 1968[). and the A-X electronic transition of CH at 

„/-v-„«„„lli nock t^u^ ,.,„..„ r< „„a r< ,.,u;„u 



432 nm jGrattoni ri985). The exceptions were C2 and C, which 
do not have bands in the UV, and for which we thus used the 
400 nm reference wavelength instead. One should also note that 
850 and 1600nm coincide respectively with the maximum and 
minimum absorbtion of the H ion, which is the most important 
contributor to the continuum opacity in red giant atmospheres in 
the optical to near-infrared wavelength range. On the other hand, 
the continuum opacity at 400 nm is dominated by the contribu- 
tion from metals. Therefore, the choice of the three wavelengths 



offers a possibility to study the interplay between the different 
sources of continuum opacity and line formation. 

The excitation potentials were selected to cover the range 
of x combinations possible for real (i.e., nonfictitious) ele- 
ments/lines. We used^ = 0-6 eV for the neutral atoms, 0- 10 eV 
for the ions (in both cases with a step of Ax = 2 eV), and -4 eV 
(A^- = 1 eV) for the molecules. Exceptions were Oi and Ni: 
their real lines are characterized by very high excitation poten- 
tials (x > 9 eV) and thus the x range for these two species was 
chosen to be identical to that of the ions. 

We finally stress that some of the elements/species stud- 
ied here have only a few real lines that can be used in the 
spectroscopic diagnostics of red giant stars (e.g., Lii, Ni, Oi, 
S 1). Moreover, we also included several elements whose spec- 
tral lines are inaccessible for observations in real life (e.g., Kn, 
S n). We stress that all these elements were included to help un- 
derstand the trends and properties of spectral line formation in 
the presence of realistically modeled convection, and to identify 
the physical causes behind them. Obviously, the exact values of 
abundance corrections for the combinations of atomic parame- 
ters where spectral lines of these elements do not exist (or may 
not be observed in stellar spectra) can only be of academical in- 
terest. 



2.2.2. Spectral line synthesis with 3D, (3D), and 1D model 
atmospheres 

Three-dimensional spectrum synthesis calculations are very 
time-consuming and thus performing them by using the entire 
3D model sequence would be impractical. To make the task man- 
ageable, synthetic spectral line profiles were computed using 14 
3D model structures (snapshots) selected from the sequence of 
seventy 3D models, fully relaxed to a quasi-stationary state (see 
Sect. 12.1. H i. To speed-up calculations further, spectral synthesis 
was carried out on a coarser x, y grid of 50 x 50 points, i.e. using 
only one third of the grid points of the original 3D model box in 
each horizontal direction. Full-resolution test calculations per- 
formed on the original 150 x 150 grid show that the reduced 
horizontal resolution has a negligible effect on the properties of 
synthesized spectral lines. 

For all elements investigated here, spectral line profiles cor- 
responding to the 3D hydrodynamical model were calculated 
for each individual 3D model structure (snapshot) in the 14- 
snapshot ensemble. A composite 3D line profile was then con- 
structed by co-adding the line profiles corresponding to all four- 
teen snapshots. Similarly, (3D) line profiles were calculated us- 
ing a sequence of (3D) models obtained according to the pre- 
scription given in Sect. 12.1.21 Spectral line profiles were calcu- 
lated for each of the fourteen (3D) models and then co-added to 
produce a composite (3D) line profile. A microturbulence veloc- 
ity of ^ m i c = 2.0km/s was used in the line synthesis calculations 
with (3D) and ID models. The choice of £ m i c is not critical since 
we used only weak unsaturated lines in the following (Sect. 12731 . 

To enable a strictly differential comparison between the pre- 
dictions of the 3D and ID models, 3D, (3D), and ID spectrum 
synthesis computations were made with the same spectrum syn- 
thesis code, Linfor3E0. Additionally, as we stated above, the 
3D and ID models used in the computations shared identical at- 
mospheric parameters, chemical composition, equation of state, 
and opacities. We therefore tried to minimize the differences 
in the model calculation procedure and line synthesis compu- 
tations, so that any discrepancy in the predictions obtained with 



http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 
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Fig. 1. Top panel: temperature profiles in the three model atmo- 
spheres of a red giant: 3D hydrodynamical (gray scales showing 
the logarithmic probability density of the 14-snapshot ensem- 
ble), average (3D) (red solid line, average over the 14-snapshot 
ensemble), and ID LHD model atmosphere with o;mlt=1-0 (red 
dashed line). Horizontal bars indicate approximate formation re- 
gions of Fei, Fen, and CH lines in the 3D and ID models, at 
different wavelengths and line excitation potentials (bars mark 
the regions where 90% of line equivalent width is acquired, i.e., 
between 5% and 95% in the cumulative line depression contribu- 
tion function, see Appendix IB! . Bottom panel: RMS horizontal 
temperature fluctuations at constant TR OSS in the 3D model (black 
solid line, 14-snapshot ensemble); and the difference between 
the temperature profiles corresponding to the average (3D) (14- 
snapshot ensemble average) and the ID model (blue dashed 
line). 

3D and ID models could be traced back to the differences in 
physical realism invoked in the two types of models. 

2.3. 3D-1D abundance corrections 

The influence of convection on the spectral line formation and 
the resulting line strengths was investigated with the aid of 3D- 
1D abundance corrections. The 3D-1D abundance correction, 
Asd-id, was defined as the difference in the abundance A(X[) of 
the element X; obtained for a given equivalent width of a par- 
ticular spectral line with the 3D hydrodynamical and classical 
ID model atmosph eres, A3D-1D = ^(Xi)3D — A(X;)id (see, e.g., 
Caffau et al.ll201 lb . The contribution to the 3D— ID abundance 
correction comes from two major constituents: (a) the correc- 
tion due to the horizontal temperature fluctuations in the 3D 
model, A3d-<3d> = ^4(Xj)3D - A(X;)(3d>, and (b) the correction 
due to differences between the temperature profiles of the av- 
erage (3D) and ID models, A <3D) _i D = A(Xi) <3D) - A(Xi)i D . 
The full abundance correction is a sum of the two constituents, 

A3D-ID = A3 D _(3D) + A(3D>-1D- 

The 3D-1D abundance corrections were always calculated 
for weak lines (equivalent width < 0.5 pm). The reason for this 
choice was that these weak lines are supposed to be on the linear 




Fig. 2. Number densities of Fe 1, Fe n, and CH (top-down), plot- 
ted versus TR OSS for the three model atmospheres of a red gi- 
ant: 3D hydrodynamical (gray scales showing the logarithmic 
probability density of the 14-snapshot ensemble), average (3D) 
(solid line, 14-snapshot ensemble average), and ID LHD model 
with Q'mlt— 1-0 (dashed line). The number densities of Fei and 
Fen are provided as fractions of the total iron number den- 
sity, n(Fe tot ), whereas that of CH is given on a scale where 
n(H tot ) = 12.0. 



part of the curve-of-growth, where their equivalent width is inde- 
pendent of the microturbulence velocity, fmic, used with the (3D) 
and ID models. Hence, the derived 3D-1D abundance correc- 
tions become independent of the choice of the microturbulence 
parameter. 

3. Results and discussion 

3.1. General properties of spectral line formation in the red 
giant atmosphere 

Convection indeed plays an important role in shaping the upper 
atmosphere of the red giant studied here: convective up-flows 
and down-drafts dominate close to the optical surface (optical 
depth tr oss ~ 1), whereas the shoc k- wave activity is most promi - 
nent in the outer atmosphere (see iLudwig & Ku cinskas 20121) . 
This alters the thermal structure of the atmosphere, in particular, 
by causing spatial and temporal variations of the temperature 
profiles in the 3D model (Fig.[T). Since local temperature sets 
the physical conditions for spectral line formation, it is reason- 
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Fig. 3. Abundance corrections for spectral lines of neutral atoms, plotted versus the difference between their ionization energy and 
line excitation potential, Ei on ~x- Three types of abundance corrections are shown: A3d-(3d> (left column), A(3d>-id (middle column), 
and A3D-1D (right column). Abundance corrections are provided at three different wavelengths: 400 nm (top row), 850 nm (middle 
row), and 1600 nm (bottom row). The ionization energies of neutral atoms used in plotting this figure are provided in TableQ] 



able to expect that line profiles predicted by the 3D, (3D), and 
ID models will also be different. 

Obviously, there are significant differences between the 
number densities of chemical species predicted by the 3D hy- 
drodynamical and classical ID models (Fig.[2j. To a large extent, 
this behavior is defined by the atomic/molecular properties of in- 
dividual species, such as ionization and dissociation potentials in 
case of atoms and molecules, respectively. This leads to different 
sensitivities of the number densities to temperature fluctuations. 
For example, species that are in the minority ionization stage at a 
given depth in the atmosphere (e.g., Fe 11 at log tr oss < -0.2) are 
very sensitive to temperature fluctuations, since small changes 
in the degree of ionization lead to a large spread in their number 
densities (as indicated by the width of the density plot in Fig. [2]). 

For a given set of spectral line parameters (such as wave- 
length, excitation potential, and oscillator strength), the line for- 
mation region is essentially defined by the temperature. Since 
differences between the temperature profiles of the average (3D) 
and ID models are small in the entire range of optical depths 
relevant to the line formation (|Ar| < 200 K in the range of 
logTR OSS < 1.0, Fig.Q]), this leads to very similar line formation 
spans in the (3D) and ID models. On the other hand, the pres- 
ence of horizontal temperature fluctuations in the line forming 
layers plays a major role in the 3D line formation process, lead- 
ing in general to line strengthening with respect to the (3D) case. 
The amplitude of the deviations from the average (3D) tem- 
perature profile, as defined by the RMS horizontal temperature 

fluctuations (AFrms = ^{(T - 7o) 2 ). tJ -,/, where {.) x ,y,t denotes 
temporal and horizontal averaging on surfaces of equal optical 
depth, and To = {T) x>y>t , is the depth-dependent average temper- 
ature), is monotonically decreasing throughout the entire photo- 



Table 1. Ionization energies of various neutral atoms. 



Element 


F eV a 


Element 


F- eV a 


Element 


F eV a 


Li I 


5.39 


SI 


10.36 


Mnl 


7.43 


Nal 


5.14 


KI 


4.34 


Fel 


7.90 


Mgl 


7.65 


Cal 


6.11 


Col 


7.88 


All 


5.99 


Til 


6.83 


Nil 


7.64 


Si I 


8.15 


CrI 


6.77 


Znl 


9.39 



a NIST database, https://www.nist.gov 



sphere, from 500 K at log tr oss = + 1 .0 to 50 K at log tr oss = -5 
(Fig-[D- We thus expect the differences between 3D and (3D) 
line formation to show up most clearly for high-excitation lines, 
forming in the deep photosphere where the horizontal temper- 
ature fluctuations are large. At the same time, the differences 
between the average (3D) and ID temperature profiles will also 
be most pronounced in this part of the atmosphere, such that the 
(3D)-1D effects should also be strongest for the high-excitation 
lines. 

3.2. Abundance corrections for lines of neutral atoms 

The abundance corrections obtained for the bulk of neutral atoms 
are plotted in Fig. [3] versus the difference between their ion- 
ization and excitation potentials, Ei on - x- This choice of ab- 
scissa was motivated by the fact that for chemical elements 
that are almost completely ionized in the line-forming regions, 
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it is the difference £i 0n - x (and °f course the thermodynam- 
ical conditions) that defines the number density of the neutral 
atoms per unit mass. The line opacity of this minority species 
is thus proportional to the comb ined Saha -Boltzmann factor, 
K[ ~ exp{+(£ , ; on ~x)/kT} (cf . iGrayl 120051 in the analysis of 
the temperature dependence of the line strength). The abun- 
dance correction curves shown in Fig. [3] must therefore fall on 
top of each other for all neutral atoms of elements that are 
strongly ionized (i.e. those with sufficiently low ionization po- 
tential, Ei on ^ 6 eV), as explained in more detail in AppendixlAl 
Indeed, this is clearly the case for the abundance corrections of 
Li i, Na i, and K i, which are in their minority ionization stage 
throughout the entire atmosphere of this particular red giant. For 
this type of atoms, the total abundance correction, A3D-1D, and 
its constituents, A3d-<3d> and A<3d>-id, are confined to the range 
of -0.1 . . . + 0.05 dex, with comparable contributions (of differ- 
ent sign) from the A3d-(3d> and A(3d>-id corrections. All correc- 
tions are more negative in the near-IR at A 1600 nm than in the 
red at A 850 nm. 

For higher ionization potential, neutral atoms gradually turn 
into majority species, and the abundance correction curves in 
Fig.[3]begin to separate. The two type of lines behave radically 
different, and this can be clearly distinguished in the A3D-1D ver- 
sus Ei on - x plot. This behavior is simply a consequence of the 
combined action of ionization and excitation, as demonstrated 
in AppendixlAl The point is that the ionization factor dominates 
over the excitation factor as long as the neutral atoms are a mi- 
nority species, and thus the line opacity decreases with increas- 
ing temperature; low-excitation lines are then most temperature- 
sensitive. The reverse is true for a neutral majority species. 
Here the excitation factor dominates over the ionization fac- 
tor, and the line opacity increases with increasing temperature, 



ki ~ exp{-x/kT] (cf. lGravll2005l Chapter 13, Case 1). In this 
situation, high-excitation lines are most temperature-sensitive. 

The largest (most negative) 3D corrections are obtained for 
Zni and S 1, reaching down to A3d-<3d> ~ -0.20 and -0.28 dex, 
respectively. This is explained by the fact that these species have 
the highest ionization potentials of all atoms shown in Fig. [3] the 
high-excitation lines of these atoms have the most temperature- 
sensitive Boltzmann factor, and at the same time form in the 
deep photosphere were the temperature fluctuations are more 
pronounced than in the higher photospheric layers were the lines 
of the minority species originate. Surprisingly, the corrections 
for all atoms of majority type are much smaller at A 1600 nm 
than at A 850 nm (see Appendix [B] for a detailed explanation). 
We note that neutral atoms with the highest ionization poten- 
tials behave as the majority ions shown in Fig.|4] where Ni 
(Zijon = 14.53 eV) and Oi (E[ on = 13.62 eV) have already been 
included. The abundance corrections for the ions (and neutral 
majority species) show the same strong wavelength dependence. 
They are discussed in the next section. 

3.3. Abundance corrections for lines of ionized atoms 

Figure|4] displays the abundance corrections for a selection of 
ionized atoms (plus N 1 and 1), showing a rather uniform de- 
pendence on excitation potential x- We first consider the ions 
that represent the majority ionization stage. Their line opacity 
(per unit mass) depends on temperature as K[ ~ exp{~x/kT}. 
Owing to its low first ionization potential, K 11 is present in its 
majority ionization stage throughout the entire atmosphere of 
this red giant (K 11 is only plotted for tracing the behavior of ele- 
ments that are nearly 100% ionized; Kn lines are inaccessible to 
observations in real red giant atmospheres). The majority neu- 
tral atoms Ni and Oi show exactly the same temperature sen- 



A. Kucinskas et al.: Three-dimensional hydrodynamical C0 5 BOLD model atmospheres of red giant stars 



7 



3D - <3D> 



<3D> - ID 



3D - ID 




-a 

§ -0.3 _ 
3 0.1 

0.0 

-0.1 

-0.2 

-0.3 [ I I ; 

012340123401234 

Z> eV 

Fig. 5. Same as in Fig.[4]but for molecular lines. The bluest wavelength (top row) corresponds to the real wavelength of molecular 
bands in this spectral range (see legend), except for C2 and CO, which do not have bands in the UV. The curve corresponding to C2 
coincides with that of CN in the lower left panel. Note the different y-scale in the lowest three panels. 




sitivity (see Sect. 13.21 ). and hence their 3D-1D abundance cor- 
rections coincide with those of the majority ions. The correc- 
tions vanish for ground-state lines and increase steadily in am- 
plitude toward the high-excitation lines. At^ = 10 eV, A3d-<3d> 
amounts to ss -0.6 dex, while A<3d)-id reaches w +0.2 dex, 
hence A3D-1D ~ -0.4 dex. For ions that are no pure majority 
species, the abundance corrections depend weakly on the ion- 
ization potential, systematically increasing in absolute size with 

^lon- 

Basically, the curves shown in Fig. [4] reflect the temperature 
sensitivity of the line opacity, which simply scales with the ex- 
citation potential, d log K(/d log T = %/(kT) (valid for the ma- 
jority species). All abundance corrections are therefore close to 
zero for x — 0, and increase systematically with x- Since the 
(3D) model is cooler in the line-forming regions than the ID 
model, the same line is weaker in the (3D) model than the ID 
model, such that A(3d>-id ~ ~x(\°?>Tqd) - log T^) is posi- 
tive. On the other hand, A3d-(3d> is negative, because the hor- 
izontally averaged line opacity in the 3D model exceeds the 
line opacity obtained for the horizontally averaged temperature, 
(k() > K(((T)) « Kt (1<3D)), essentially because of the highly non- 
linear temperature dependence of K{ (d 2 K(/dT 2 > 0). 

In this simple picture, the abundance corrections should be 
similar for all wavelengths, since the temperature sensitivity of 
the line opacity is wavelength independent (apart from a weak 
dependence through the stimulated emission factor). Obviously, 
this is not the case: the abundance corrections are much larger in 
the red at A 850 nm than in the near-IR at A 1600 nm, even though 
the near-IR lines form at deeper photospheric layers where 
the amplitude of the temperature fluctuations should be larger 
(Fig-U lower panel). A detailed explanation of this counter- 
intuitive behavior is given in Appendix FBI In short, it is not the 



line opacity alone that determines the strength of a spectral line, 
but rather the ratio of line-to-continuum opacity. It turns out that 
the continuum opacity at A 850 nm (due to H bound-free ab- 
sorption) is almost independent of temperature, such that the 
T-dependence of the line opacity dominates the abundance cor- 
rections. At A 1600 nm, however, the continuum opacity (due to 
H free-free absorption) is strongly T-dependent, leading to a 
substantial reduction of the temperature sensitivity of the ratio 
kc/k c , and hence to much smaller abundance corrections than at 
/I850nm. 

Finally, we note that the S 11 line exhibits by far the largest 
abundance corrections, reaching A3d-<3d) ~ —1.1 dex at x - 
10 eV. This outstanding behavior is related to the high ion- 
ization potential of Si (.Eion = 10.36 eV), making Sua true 
minority ion. In this ca se, the line opacity is proportional to 
exp{-(£ ion +x)/kT} (cf. lGravll2005l Chapter 13, Case 3). This 
extremely high temperature dependence gives rise to the extraor- 
dinary abundance corrections. In compliance with this interpre- 
tation, A3d-<3d> for S 11 atx = is already as large as A3d-<3d> at 
X — 10 eV for the normal majority ions. We note, however, that 
this is probably a purely academic case, since such lines would 
certainly be too weak to be observable. 



3.4. Abundance corrections for molecular lines 

Figure|5] displays the abundance corrections for a selection of 
molecules as a function of the excitation potential x- The to- 
tal abundance corrections (right column) range from A3D-1D ~ 
-0. 16 dex for C 2 at A 850 nm to « +0. 1 dex for NH and OH in the 
UV. As in the case of neutral atoms, the A3d-<3d> an d A<3d>-id 
corrections are often of opposite sign, making the total abun- 
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dance correction, A3D-1D, slightly smaller than the granulation 
correction A3d-<3d>- 

A detailed interpretation of the molecular abundance cor- 
rections is more complicated than in the case of neutral atoms 
and ions, because the dissociation equilibria of all molecules are 
closely coupled, such that their number densities are intimately 
related and in some cases show an unexpected temperature de- 
pendence. For basic orientation, Fig. IC.ll (Appendix [Cb shows 
the depth-dependence of the molecule concentrations (number 
densities normalized to total number density of carbon nuclei) 
in the ID LHD model. 

First of all, because of its high dissociation energy of Do = 
11.1 eV, CO is a majority species throughout the whole photo- 
sphere; practically all carbon nuclei are locked in this molecule 
at tr oss < 1, whereas CO dissociates quickly in the deeper pho- 
tosphere. In terms of abundance corrections, CO behaves like the 
majority species shown in Fig. [4] The dependence on excitation 
potential is weaker, however. Owing to dissociation in the deeper 
layers, the excited CO lines form higher up in the atmosphere 
than the majority ions with the same^-, and thus feel weaker tem- 
perature fluctuations and a reduced (3D)-1D temperature differ- 
ence. As a consequence, the 3D abundance corrections are small 
in the studied range of excitation potentials. 

The molecules NH and OH belong to a different class of 
molecules. Their constituents (H, N, and O) are majority species 
but only a tiny fraction of N and O is bound in NH and OH 
molecules, respectively. In this situation, the molecular line 
opacity is proportional to the Boltzmann factor exp{(Do-x)/kT], 
i.e. it depends on the difference between dissociation energy Do 
and excitation potential^. The largest abundance corrections are 
therefore found for x — 0, where the line opacity has the high- 
est temperature sensitivity (see Appendix ICl Eq. IC.5t . In fact, 
the abundance corrections of these molecular lines behave very 
much like those of the true minority neutral atoms shown in 
Fig,[3] This is also true for their wavelength dependence. 

The molecules CN, CH, and C2 fall in yet another category. 
Their temperature dependence is opposite to that of the other 
molecules, i.e. their number density increases toward higher 
temperature (at constant density or constant pressure). The rea- 
son is a strong coupling to the dissociation equilibrium of CO, 
which controls the number density of free carbon atoms. Higher 
temperatures lead to some dissociation of CO, hence an in- 
creased concentration of carbon atoms, and in turn to higher den- 
sities of CN, CH, and C2. CH is more sensitive to this coupling 
than CN, because the latter has a higher dissociation energy. The 
C2 molecule is most sensitive since its concentration depends 
on the square of the number density of free carbon atoms. The 
abundance corrections A3d-<3d> at A 850 nm are therefore largest 
(most negative) for C2. 

At A 1600 nm, the fluctuation of the continuum opacity are 
stronger than in the red at A 850 nm, because (i) the tempera- 
ture fluctuations are larger in the deeper layers where the near- 
IR lines form, and (ii) the temperature dependence^ of the con- 
tinuum opacity is stronger in the near-IR (see Fig. IB. 21 ). For 
molecules like NH and OH, the 3D-(3D) corrections are larger 
(more negative) at A 1600 nm since K( decreases while k c in- 
creases with T, such that the fluctuations of the ratio K(/k c are 
stronger than the fluctuations of K[ . The molecules CN, CH, and 
C2, on the other hand, have the opposite temperature dependence 
of K(, and thus the fluctuations of the continuum opacity effec- 
tively cancel the fluctuations of the line opacity. The corrections 



A3d-<3d> for CN, CH, and C2 are therefore much smaller in the 
near-IR than in the red. 



3.5. Abundance corrections at 2.5 fim 

As a test case, we also calculated the 3D-1D corrections for 
several fictitious lines of different neutral, ionized atoms and 
molecules at 2.5 fim (not shown in Figs. [3}|5j. In this case the 
A3d-<3d>, A<3d)-id, and A3D-1D abundance corrections, as well 
as their dependence on the line parameters, were very similar to 
those obtained at 1 .6 fim. This result is hardly surprising since 
opacity in the red giant atmospheres at 1.6 - 2.6 fim is domi- 
nated by H free-free transitions, the efficiency of which varies 
very little in this wavelength range. Thus, the line formation re- 
gions are essentially the same at both wavelengths, which leads 
to very similar line formation properties and 3D-1D abundance 
corrections. 



3.6. Line formation in the 1D model atmospheres: influence 
of the mixing-length parameter ff MLT 

According to the classical Schwarzschild criterion, ID LHD 
models of the red giant studied here are convectively stable 
in the outer photospheric layers (logT R OSS < 1), irrespective 
of the mixing length param eter, omlt JKucinskas et al. 2005; 



including the indirect temperature dependence via electron pres- 



iLudwig & Kucinskasll20l2h . Since line formation depths typi- 
cally never reach deeper than log tr oss =s 0.6, the choice of ^mlt 
should have no effect on the resulting spectral line strengths. The 
only exceptions may occur in case of near-infrared lines and/or 
lines characterized by very high excitation potentials (> 8 eV), 
i.e. those that form deepest in the photosphere and thus may 
be sensitive to differences in the temperature profiles computed 
with different mixing-length parameters. 

To verify whether this can actually happen, we ran several 
test calculations using fictitious Fe 1 (0, 4 eV), Fe 11 (0, 4, 8 eV), 
and On (0, lOeV) lines located at 850nm and 1600nm. We 
found that differences between the abundances derived using 
model atmospheres calculated with different mixing-length pa- 
rameters (ffMLT = 1,2) were always very small. In fact, they 
never exceeded 0.02 dex and thus may be safely ignored in most 
situations related with the abundance work in this particular red 
giant. 

3. 7. How to use the 3D abundance corrections 

We recall that in this work the 3D-1D corrections were calcu- 
lated for very weak spectral lines. While such corrections pro- 
vide a convenient way to estimate the importance of 3D hydro- 
dynamical effects in the spectral line formation, abundance cor- 
rections for stronger lines may be different because of saturation 
effects, leading to differences in formation depths and introduc- 
ing a sensitivity to velocity fields (see Appendix lB.3t . 

The theoretical 3D abundance corrections were derived from 
the comparison of the equivalent width of the same (artificial) 
spectral line computed with a 3D model and a ID model. Since 
the 3D and the ID models use the same stellar parameters, 
atomic data, and numerical methods (as far as possible), the re- 
sulting differential 3D corrections should be applicable to any 
ID model, irrespective of the physical details. 

For given stellar parameters (T e ff, logg, [M/H]) and a given 
element and ionization stage, the 3D abundance correction de- 
pends only on the energy of the lower level of the transition, the 
wavelength of the line, and in principle also on the mixing-length 
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parameter used for the ID model, A3D-iD0f , -i> »mlt), provided 
that the line is weak (on the linear part of the curve-of-growth). 
In this case, the corrections given in Figs. [3] @] and [5] can be 
readily applied to the ID LTE abundance determinations per- 
formed with any standard ID mixing-length model atmosphere: 
■A(Xi)3D = A(X;)id + A3d-id- As shown in Sect. 13.61 the depen- 
dence on q-mlt is negligible for the present stellar parameters. 

In the general case of stronger lines, the 3D corrections de- 
pend in addition on the the equivalent width of the line, W, and 
on the microturbulence parameter used for the ID model, 
A3D-iD0f, A, £mi C ; (Emlt), as illustrated in Fig. IB. 41 for Fen, 
X — 10 eV. Except for this example, we do not provide the de- 
pendence of A3D-1D on W and £ m j c in this work. The given cor- 
rections for the weak line limit can then only serve as an indica- 
tion of the line's susceptibility to 3D effects. Figure lB~4"l suggests 
that the 3D correction for stronger lines is always more positive 
than in the weak line limit. The figure also seems to indicate that 
A3D-1D becomes independent of W for strong lines, suggesting 
that a 'strong line limit' might be a useful quantity. Additional 
investigations are needed to see whether these properties remain 
valid in general. 

4. Conclusions 

We have investigated the influence of thermal convection on the 
spectral line formation in the atmosphere of a red giant located 
close to the RGB tip (T eS = 3660, log g = 1 .0, [M/H] = 0.0). To 
this end, we synthesized a large number of fictitious atomic and 
molecular lines of astrophysically important tracer elements, us- 
ing for this purpose the 3D hydrodynamical C0 5 B0LD and clas- 
sical ID LHD stellar model atmospheres (both types of models 
sharing identical atmospheric parameters, opacities, equation of 
state, and chemical composition). The influence of convection 
on the line formation was investigated by focusing on the dif- 
ferences between the abundances inferred from a given (weak) 
spectral line by the requirement of producing a given equivalent 
width with the 3D and ID model atmospheres, assuming LTE. 

Overall, convection plays a considerable role in the atmo- 
sphere of the red giant studied here. There are significant hori- 
zontal temperature fluctuations seen at different optical depths in 
the atmosphere that are caused either by convective up- or down- 
flows (inner atmosphere), convective overshoot, and/or shock 
wave activity (outer atmosphere). This leads to substantial hor- 
izontal variations in the number densities of certain chemical 
species, causing significant deviations from the predictions of 
a classical ID model. 

Spectral line formation is more or less strongly affected by 
convection in this particular red giant, depending on the line pa- 
rameters and the chemical species under consideration. The dif- 
ferences in elemental abundances inferred from a given spectral 
line with the 3D and ID model atmospheres, A3D-1D, are most 
pronounced for high-excitation lines of ions and atoms of pre- 
dominantly neutral elements. Their abundance corrections grow 
larger with increasing excitation potential, reaching values of 
A3D-1D ~ _ 0.4 dex for excitation potential x = 10 eV. The 
main physical reason for this 3D-1D difference is the increas- 
ingly nonlinear temperature dependence of the line opacity as x 
increases. In a 3D atmosphere, the horizontal temperature fluc- 
tuations then lead to an enhancement of the effective line opacity 
with respect to the ID case. 

Lines of neutral atoms of predominantly ionized elements 
show significantly smaller corrections, with A3D-1D not exceed- 
ing +0.1 dex. Here the temperature dependence of the line opac- 
ity due to ionization and excitation tend to cancel partially. 



Molecular lines show a more complex behavior due to the 
strong coupling between the different dissociation equilibria. CO 
is least susceptible to 3D effects, with |A3d-idI < 0.03 dex. 
NH and OH represent a different category of molecules; their 
abundance corrections fall in the range |A3d-idI < 0.1 dex. 
The largest corrections are found for C2, reaching as low as 
A3D-1D ~ -0.16 dex. 

It is important to emphasize that the 3D-1D corrections for 
all chemical species studied here show a significant wavelength 
dependence. In most cases, the abundance corrections are signif- 
icantly smaller in the near-infrared at A 1600 nm than in the opti- 
cal spectral range. Careful investigation reveals that this does not 
indicate that the atmospheric layers where the infrared lines orig- 
inate are least affected by convection. Rather, the strong wave- 
length dependence of the 3D corrections is related to the fact 
that the continuum opacity is much more temperature depen- 
dent around A 1600 nm (mainly H free-free absorption) than 
around A 850 nm (mainly H bound-free absorption). It should 
be stressed though that this conclusion was reached for the par- 
ticular red giant studied here and may change with the atmo- 
spheric parameters of the star under investigation. 

One may thus conclude that the spectral line formation in 
a red giant photosphere is a delicate process, governed by the 
subtle interplay between microscopic (atomic line parameters) 
and macroscopic (local temperature fluctuations) physics. The 
size of the 3D-hydrodynamical fluctuations and their effects on 
the line formation process can only be assessed by dedicated ra- 
diation hydrodynamics simulations and detailed 3D line forma- 
tion calculations. For the present example, we find substantial 
3D abundance corrections, especially for ionized atoms, sug- 
gesting that even at solar metallicity, the accuracy of ID LTE 
stellar abundance work can be significantly improved by exploit- 
ing the additional information provided by realistic 3D model 
atmospheres. 
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Fig.A.l. Abundance correction A3d-<3d>, computed according 
to Eqs. ( IA.U - dA.41 . versus the difference between ionization 
and excitation potential, E\ ari - x- Each curve corresponds to a 
different value of E\ ori . The thermodynamic variables 9{x,y) and 
p e (x, y) were taken from the 3D model at monochromatic optical 
depth log rg5o = -0.48, where the mean temperature is 3360 + 
140 K. The temperature dependence of the partition functions 
Uq and U\ has been neglected. 
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Appendix A: Analysis of the abundance corrections 
for lines of neutral atoms 

The temperature sensitivity of spectral lines originating from the 
neutral atoms of partially ionized species (e.g. Mg i, Ca i, Fe i) is 
governed by the Saha and Boltzmann equations, i.e. by changes 
of the degree of ionization and of the excitation of the line's 
lower level. In the following, we develop a simplified model of 
the 3D-(3D) abundance corrections, which result from horizon- 
tal fluctuations of the thermodynamical quantities. 

The ratio of the total number of neutral atoms no (with ion- 
ization energy E mn and partition function Uq) to the total number 
of singly ionized atoms n\ (with partition function U\) at tem- 
perature T = 5040/ and elec tron pressure p e is given by Saha's 
equation (see e.g. Gray 2005, Eq. 1.20) 



, «o 
log — 

«i 



E ion 8 + 2.5 log 6 + log p e + log - 0.08 



The opacity of a line transition with lower level excitation po- 
tential x is proportional to the number of absorbing atoms (per 
unit mass) in this state of excitation, 



no 



1 



l+t\0, Pe ) 



(A.2) 



In the presence of horizontal fluctuations of 6 and p e , the average 
line opacity is amplified with respect to the (3D) case by a factor 



(K((6, p e ))x,y 
K{ ((0)x,y, (Pe)x,y) 



(A3) 



where {.) x ,y denotes horizontal averaging at constant continuum 
optical depth t c . Assuming that fluctuations of the continuum 
opacity and the source function can be neglected, the abundance 
correction for weak lines can be estimated as 



^3D-<3D> 



(A.4) 



Figure IATTI shows the result of computing A3d-<3d> according to 
Eqs. dA.U - dA.4t for different combinations of^ and E mn . The 
figure illustrates how the curves A3d-<3d> versus E ion ~x change 
systematically as the parameter E mn increases from to 9 eV. For 
Ei 0B ^ 6 eV, the neutral atoms are a minority species, / « 1, and 
we see from Eq. ( IA.2t that in this case 



f 10~ 9x ~ io e(£ion ^ ) 



(A.5) 



depends only on the difference (Ei on - x), such that all curves 
with Ei on ~ 6 eV fall on top of each other. As Ei on increases 
further, the vertex of the curves moves from (Ei on ~x)max = eV 
to S3 4 eV to 9 eV as the ionization balance shifts from {f) x ,y ~ 
(Eion = eV) to * 1 (E^n = 7 eV) to > 1000 (E ion = 9 eV). In 
fact, it can be shown analytically that 



(-Eion ~X)n 



</>* 



1 + </>*,, 



(A.6) 



= log f(0,p e ). 



(A.l) 



Admittedly, the description of the abundance corrections de- 
veloped above is severely simplified. It ignores the fact that the 
line formation region is extended and that the location of its cen- 
ter of gravity depends sensitively on Ei OB and x (cf. Figs. IA.2l 
IB. II IC.2l i. Also, fluctuations of the continuum opacity and the 
source function were neglected. Nevertheless, the systematics 
seen in Fig. I A. II provides a basic explanation of the detailed nu- 
merical results presented in Sect. 13. 21 Fig. [3] (especially middle 
left panel). 

Appendix B: Analysis of the abundance corrections 
for high-excitation lines of ions 

In the following, we analyze in some detail the abundance cor- 
rections derived for the high-excitation Fe 11 lines, which are rep- 
resentative of the ionized atoms and show the largest 3D correc- 
tions (see Fig.|4]l. Evaluating the equivalent width contribution 
functions of this line in the 3D and the ID models, we determine 
the physical cause of the abundance corrections. In particular, 
we can understand the sign of the 3D-(3D) and (3D)-1D cor- 
rections and explain why these corrections are so much smaller 
at A 1600 nm than at A 850 nm. 

For simplicity, we consider only vertical rays (disk-center 
intensity) in a single snapshot from the 3D simulation, noting 
that the qualitative behavior of the abundance corrections is sim- 
ilar for intensity and flux, and does not vary much in time, i.e. 
the 3D-(3D) correction is always strongly negative, while the 
(3D)-1D correction is slightly positive at A 850 nm. 
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Fig.A.2. Disk-center (yu = 1) equivalent width contribution 
functions, S(logr c ), of a weak (artificial) Fei line with excita- 
tion potential x — eV (top) and y — 5 e V (bottom), at wave- 
lengths A 850 nm, evaluated according to the weak line approx- 
imation (Eqs. IB.5llB.7l i. for a single snapshot of the 3D model, 
the corresponding (3D) average model, and the associated ID 
LHD model used in this work. The contribution functions, orig- 
inally defined on the monochromatic optical depth scale, have 
been transformed to the Rosseland optical depth scale. 

B.1. Formalism 



Following Magain ( 1986), the line depression contribution func- 
tion (for vertical rays, LTE), Cf , is defined in Linfor3D as 

Cf (AA, T c ) = ( T](AA, t c ) u c (t c ) exp{- (r c + t ( (AA, t c ))} ) xy 

(B.l) 

where t c is the continuum optical depth, rj{AA, r c ) = 
K((AA, t c )/k c (t c ) is the ratio of line opacity to continuum opac- 
ity, m c (t c ) = 7^(t c ) - S c (t c ) is the difference between outgoing 
continuum intensity and source function^, and (t c + T((AA)) is 
the total optical depth in the line; angle brackets (.) A y indicate 
horizontal averaging at constant continuum optical depth. Note 
that rj (and T( ) vary with wavelength position in the line profile, 
AA, whereas r c and u c can be considered as constant across the 
line profile. Then the absolute line depression at any wavelength 
in the line profile is 



D/(r c = 0, AA) 



f 

Jo 



Cf (AA, r' c )dT' c 



(B.2) 



Cf(r c ) 



Cf (AA', 

oo 



Defining further the equivalent width contribution function as 

t c ) d AA' , (B.3) 
the equivalent width of the line is finally computed as 

i r°° 

Vc) JO 

= TTT f S(log<)dlog/ c , (B.4) 

where (I c ) is the horizontally averaged emergent continuum in- 
tensity, and S is defined as 



S(t c ) = lnl0T c Cf(T c ). 



(B.5) 



In the limit of weak lines, we can assume that T( <zz t c over 
the whole line formation region, and Eqs. ( IB. II ) and ( IB- 31 > sim- 
plify to 



Cf (AA, t c ) = exp{-T c } ( ?7(AA, t c ) k c (t c ) ) t 
and 

Cf(T c ) = exp{-T c } < 77o(T c ) m c (t c ) ) c ,, , 
where 



%(t c ) 



f 



j](AA', r c ) d AA' . 



(B.6) 
(B.7) 

(B.8) 



4 u c is also called 'source function gradient' since u c = dS c /dr c in 
the diffusion approximation, and if 5 c is a linear function of t c . 



In this weak line limit, Cf , and hence Cj are strictly propor- 
tional to the line opacity, and the equivalent width scales lin- 
early with the g/-value of the line (or the respective chemical 
abundance). Then the 3D abundance corrections can simply be 
obtained from the equivalent widths as 

A 3 d-id = - log {Wml Wm) , 

A 3 D-<3D> = - log (W3D/W<3D>) . 

A <3D )-id = -log(W (3D) /W 1D ). (B.9) 

B.2. Analysis of mixed contribution functions 

Figure IB. II compares the equivalent width contribution func- 
tions, S, of a weak (artificial) Fe n line with excitation potential 
X — 10 eV at two wavelengths, A 850 nm (top) and A 1600 nm 
(bottom), for the 3D model, the corresponding (3D) average 
model, and the associated ID LHD model used in this work. 
For each of the different models the area below the correspond- 
ing curve is proportional to the equivalent width of the emerg- 
ing line profile. At A 850 nm, the equivalent width produced by 
the 3D model is significantly larger than that of the ID LHD 
model, which in turn is significantly lager than that of the (3D) 
model, assuming the same iron abundance in all cases. The abun- 
dance corrections derived with Eq. ( IB. 9b are A 3 d id = -0.33, 
A 3 d-< 3 d) = -0.46, and A( 3 d>-id = +0.13 dex. These num- 
bers are fully consistent with the results shown in Fig. [4] At 
A 1600 nm, on the other hand, the equivalent widths obtained 
from the three different model atmospheres are obviously very 
similar; all abundance corrections are much smaller than those 
at A 850 nm, again in basic agreement with the results shown in 

Fig.a 

Figure lBTI also shows that the formation region of this high- 
excitation Fe ii line is well confined to a narrow region in the 
deep photosphere, mainly below tr oss = 1, where the degree of 
ionization of iron is changing rapidly with depth. As expected, 
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Table B.l. Abundance corrections for the Fen line (x — 10 eV) 
derived from mixed contribution functions &ijjc at A 850 and 
1600 nm. 



Case 


i jkm 






i jkm 




AiJ,i,m 






850 nm 


1600 nm 




850 nm 


1600 nm 





2222 


0.000 


0.000 


1111 


0.000 


0.000 


1 


3222 


0.014 


0.036 


2 111 


0.062 


-0.053 


2 


2322 


-0.451 


-0.146 


1211 


0.140 


0.150 


3 


3322 


-0.558 


-0.165 


221 1 


0.202 


0.093 


4 


2232 


0.003 


-0.132 


1121 


-0.077 


-0.130 


5 


3232 


0.046 


-0.040 


2121 


-0.015 


-0.184 


6 


2332 


-0.374 


-0.015 


1221 


0.069 


0.026 


7 


3332 


-0.460 


-0.009 


2221 


0.130 


-0.031 
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Fig. B.l. Disk-center (// = 1) equivalent width contribution 
functions, S(logr c ), of a weak (artificial) Fen line with exci- 
tation potential ^ = 10 eV, at wavelengths A 850 nm (top) and 
A 1600 nm (bottom), evaluated according to the weak line ap- 
proximation (Eqs. IB.5l IB.7l i. for a single snapshot of the 3D 
model, the corresponding (3D) average model, and the associ- 
ated ID LHD model used in this work. The contribution func- 
tions, originally defined on the monochromatic optical depth 
scale, have been transformed to the Rosseland optical depth 
scale as a common reference. 



the line originates from somewhat deeper layers at A 1600 nm 
(minimum of H opacity) than at A 850 nm (maximum of H 
opacity). According to Fig.[TJ both the amplitude of the horizon- 
tal temperature fluctuations and the ID - (3D) temperature differ- 
ence increase with depth in the range - 1 < log tr oss < + 1 , such 
that they are larger in the line formation region at A 1600 nm. 
Naively, one would thus expect the amplitude of the abundance 
corrections A3d-<3d> and A(3d>-id to be larger at A 1600 nm than 
at A 850 nm. However, this reasoning obviously fails. As we have 
seen before, the abundance corrections are found to be strikingly 
smaller at A 1600 nm. Additional analysis is necessary to resolve 
this apparent contradiction. 

To understand the origin of the abundance corrections, we 
need to understand the role of the different factors that make up 
the contribution function S, essentially u c and rjo = K(/k c . The 
different behavior of these factors in the different types of models 
determines the sign and amplitude of the abundance corrections. 
In the following, subscripts 1, 2, and 3 refer to the ID LHD 



model, the (3D) model, and the 3D model, respectively. With 
this notation in mind, we define the mixed contribution functions 



hj,k(T c ) = In 10 t c exp{-r c } ( w Ci j(T c ) 



K(J (T c ) 
Kc,k (jc) 



(B.10) 



where i = 1 . . . 3, j = 1 . . . 3, k = 1 ... 3. The three mixed contri- 
bution functions with three identical subscripts i = j = k are 
thus the normal contribution functions for the ID, (3D), and 
3D model, respectively. From each of the Syjc we can compute 
an equivalent width according to Eq. ( IB.4t . which we denote as 
Wij,k. The equivalent widths can then be used to derive abun- 
dance corrections via 



Tog(W/,;,*/W m , m , m ) . 



(B.ll) 



The numerical evaluation of the relevant abundance corrections 
is compiled in Table IB. 2l 

B.2.1 . (3D)-1 D abundance corrections 

With the help of Table IB. 21 the physical interpretation of the 
A<3d>-id abundance correction is straightforward. Columns (6) 
and (7) show the effect of the different factors that contribute to 
the (3D)-1D abundance correction. Owing to the different ther- 
mal structure of the two model atmospheres, all three factors, 
u c , Ke, and k c change simultaneously, and the full correction is 
A<3d>-id = ^2,2,2,1, listed in the last row of the Table as case (7). 
The other cases (l)-(6) refer to 'experiments' where only one 
or two of the factors are allowed to change while the remaining 
factors are fixed to expose the abundance corrections due to the 
individual factors. Case (1), for example, shows the correction 
that would result for fixed opacities, K((t c ), k c (t c ), accounting 
only for the differences in the source function gradient u c . Case 
(6) shows the complementary experiment where u c is fixed and 
both opacities are changing in accordance with the different ther- 
modynamical conditions. 

At /1850 nm, the continuum opacity is dominated by H 
bound-free absorption, which shows its maximum at this wave- 
length. The high-excitation Fe n line forms around log t^so ~ 
+0.7 (logTR OSS » +0.6), i.e. significantly below continuum opti- 
cal depth unity. At this depth, both the temperature and the tem- 
perature gradient are slightly lower in the (3D) model than in the 
ID model. As a consequence, both u c and k c decrease toward the 
(3D) model, approximately by the same factor (see cases 1 and 
4), and hence their effects cancel out. The highly temperature- 
dependent line opacity (dlogKe/86 ~ -10; = 5040/T) is thus 
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Fig. B.2. Continuous opacity due to H bf at A 850 nm (solid line, 
cloud of blue dots), and due to H ff at A 1600 nm (dashed line, 
cloud of red dots), as a function of temperature in the 3D model 
(dots) and in the (3D) model (lines), respectively. The opa cities 
have b een computed accor ding to Eqs. (8.12) and (8.13) bv lGravl 
(2005); they are given in units of cm 2 per neutral hydrogen atom, 
and have been normalized to 1 at T = 6000 K. Both opacities are 
proportional to the electron pressure, and their ratio is a simple 
monotonic function of temperature. Note that K(H bi ) is almost 
temperature insensitive between T = 4000 and 5000 K. 



the dominating factor and determines the total A(3d>-id abun- 
dance correction (compare cases 2 and 7). 

At A 1600 nm, the situation is different. Here the continuum 
opacity is mainly due to H free-free absorption. The important 
difference is that the temperature sensitivity of the H free-free 
opacity is significantly higher than that of the H bound-free 
absorption (see Fig. lB.2l >. In the line formation region around 
logri6oo ~ +0.15 (logTR OSS « +0.85), the temperature sensitiv- 
ity of K( and k c is now comparable, such that the ratio of both 
opacities is nearly the same in the two models. The corrections 
due to Kt and k c are almost equal and of opposite sign (cases 2 
and 4), and hence cancel out. At the same time, the source func- 
tion gradient is very similar in both models, and thus the correc- 
tion due to m c is small (case 1). The total A(3d>-id abundance cor- 
rection is therefore significantly smaller than at A 850 nm (case 
7). 

B.2.2. 3D-(3D) abundance corrections 

The physical interpretation of the A3d-<3d> abundance correction 
proceeds along similar lines. Columns (3) and (4) of Table lB.2l 
show the influence of the different factors that contribute to the 
3D-(3D) 'granulation correction'. The three factors, w c , Kt, and 
k c differ between the 3D and the (3D) model due to the presence 
of horizontal fluctuations of the thermodynamical conditions at 
constant optical depth t c , which then lead to more or less non- 
linear fluctuations of the factors that make up the contribution 
function. The full correction, A3d-<3d> = A3 3 3 2, allows for fluc- 
tuations in all three factors and is listed in the last row of the table 
as case (7). The other cases (l)-(6) refer to 'experiments' where 
the fluctuations are artificially suppressed for one or two of the 
factors to study the impact of the fluctuations of the individual 
factors on the resulting the abundance correction. For example, 
case (5) shows the correction that would result if the fluctuations 
of the line opacity, K({t c ), were suppressed, case (2) shows the 
complementary experiment where only K((t c ) is allowed to fluc- 
tuate, while u c and k c are fixed. 



At A 850 nm, our weak Fe n line is strongly enhanced in the 
3D model due to the nonlinear fluctuations of the line opac- 
ity. The fluctuations lead to a line enhancement, and hence to 
negative 3D-(3D) abundance corrections, whenever (iQ(T)) Xt y > 
K(({T) xy ), which happens to be the case as K[ oc exp{-E /kT) 
(roughly speaking because d 2 Kg/dT 2 > 0). As can be deduced 
from the comparison of cases (2) and (7), suppression of the 
fluctuations of both u c and k c does not change the resulting 3D 
abundance correction. We can furthermore see that the fluctu- 
ations of m c enhance the nonlinearity of the fluctuations of Kt 
(case 3) and that the fluctuations of k c diminish the nonlinearity 
of the fluctuations of K( (case 6). We conclude that the fluctua- 
tions of m c and k c must be substantial, but essentially linear, such 
that they do not produce any significant abundance corrections 
on their own (cases 1, 4, and 5). 

At A 1600 nm, the continuum opacity k c is lower than at 
A 850 nm, and our weak Fe n line forms at somewhat deeper lay- 
ers where the temperature is higher. Equally important, the tem- 
perature sensitivity of k c is distinctly higher at A 1600 nm than 
at A 850 nm, as is demonstrated in Fig. IB. 21 This fact is the key 
to understanding the drastically smaller abundance corrections 
found at A 1600 nm. 

Comparing the effect of the line opacity fluctuations for the 
two wavelengths (case 2), we see that the corresponding abun- 
dance correction is significantly smaller at A 1600 nm. This re- 
sult is unexpected, because according to Fig. Q] the temperature 
fluctuations, 6T ims , ought to be larger in the deeper layers where 
the near-IR line forms, which in turn should lead to more nonlin- 
ear fluctuations of the line opacity and hence larger abundance 
corrections at A 1600 nm compared to A 850 nm. 

Further investigations revealed that the opposite is true. The 
point is that we have to distinguish between fluctuations at con- 
stant Rosseland optical depth, tr oss , and fluctuations at con- 
stant monochromatic optical depth, t c , which are relevant in 
the present context. In fact, the higher temperature sensitivity 
of the continuum opacity at A 1600 nm reduces the amplitude 
of the temperature fluctuations at constant continuum optical 
depth T1600 with respect to the fluctuations at constant t^o, as 
illustrated in Fig. [jO] (top panel). The degree of nonlinearity of 
the line opacity fluctuations, as measured by the ratio of av- 
erage line opacity to line opacity at mean temperature, JKt = 
(k({T)) I K(({ T)), is shown in the bottom panel of Fig. IB. 31 Over 
the whole depth range, the nonlinearity of the Kg fluctuations is 
higher at constant Tg5o than at constant t^oo- Remarkably, J\( 
increases toward lower temperatures, even though the amplitude 
of the temperature fluctuations decreases with height. This is be- 
cause the temperature sensitivity of K( increases strongly as Fe 11 
becomes a minority species at lower T (cf. Fig.|2j». The fact that 
J\c is significantly higher for the red line at A 850 nm than for the 
near-IR line at A 1600 nm explains the wavelength dependence 
of the abundance corrections found for case (2), columns (3) and 
(4). 

Comparing cases (2) and (6) for the near-IR line, we see that 
the abundance correction essentially vanishes when combining 
the fluctuations of the line opacity with the fluctuations of the 
continuum opacity. We note that the fluctuations of k c are signif- 
icantly larger at log 7"i6oo = 0.15, where the near-IR line forms, 
than at log Tgso = 0.7, where the red line forms, even though the 
temperature fluctuations are lower (see Fig. IB. 3b . This is again a 
consequence of the enhanced temperature sensitivity of the con- 
tinuum opacity at A 1600 nm (Fig. lB.2l >. It thus happens that the 
abundance corrections due to the fluctuations of the continuum 
opacity and the line opacity, respectively, are comparable (com- 
pare cases 2 and 4). The net result is a cancelation of the two 
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Fig. B.3. Amplitude of the (relative) temperature fluctuations, 
STnas/iT) (top), and ratio of average line opacity to line opac- 
ity at mean temperature, 3{( = (k{(T))/k{(( T)) of a weak (ar- 
tificial) Fen line with excitation potential^- = 10 eV (bottom) 
as a function of the mean temperature (T), where angle brackets 
indicate horizontal averaging at constant monochromatic optical 
depth. Solid and dashed curves show the results of averaging on 
surfaces of constant t^q and constant T1600, respectively. Filled 
dots indicate the location of the center of the line formation re- 
gions at A 8500 and 1600 nm. 



Fig. B.4. Total 3D abundance correction A3D-1D for the artificial 
Fe 11 line with excitation potential x = 10 eV at A 850 nm (top) 
and A 1600 nm (bottom) as a function of the equivalent width 
obtained from the 3D model. The fainter (green) curves and the 
thicker (black) curves refer to the intensity (/j. = 1) and flux spec- 
trum, respectively. The abundance corrections have been com- 
puted for three different values of the microturbulence parame- 
ter used with the ID model, £ m i c = 0.0 (dotted), 1.0 (solid), and 
2.0 km/s (dashed lines). The weak line limit coincides with the 
horizontal part of the curves at low log Wid- 



effects. The total A3d-(3d> abundance correction A 1600 nm is 
therefore small. 

B.3. Saturation effects 

So far we have considered the abundance corrections for the lim- 
iting case of weak, unsaturated lines. In this limit, the abundance 
corrections are independent of the equivalent width of the line 
and of the microturbulence parameter £ m ; c chosen for the spec- 
trum synthesis with the ID models. Figure lB~4l shows how the 
results change if saturation effects are fully taken into account, 
again for the example of the high-excitation Fe 11 line. 

Obviously, the total 3D abundance correction, A3D-1D, de- 
pends strongly on both the assumed value of £ m ; c and on the line 
strength, Wjd- We notice that this holds even for very weak lines, 
and conclude that even the weakest lines used for this study are 
already partly saturated. Plotting log W3D versus log gf reveals 
that the curve-of-growth is linear, and thus saturation effects can 
be safely ignored as long as the equivalent width of the line is be- 
low W; D « 0.001 pm at A 850 nm (W; d * 0.1 pm at A 1600 nm). 
As soon as this line becomes detectable, it is no longer on the 



linear part of the curve-of-growth. This extreme behavior is of 
course related to the extreme temperature sensitivity of this high- 
excitation line, which changes the line-to-continuum opacity ra- 
tio from 77 <k 1 to rj » 1 within the line formation region. This 
is especially true at A 850 nm, where the continuum opacity is 
less temperature dependent (see above). The partial saturation 
of weak lines is not a particular property of the 3D model, but is 
seen in ID models, too. 

The top panel of Fig. lB.5l shows the equivalent width contri- 
bution functions S(tr oss ) of the same weak Fe n line (x — 10 eV, 
A 850 nm) as in Fig. IB.ll (top), but now including the satura- 
tion factor exp{-Tc} (see Eq. lB.U . Comparison with Fig. IB. II 
(top) demonstrates that including the saturation factor reduces 
the equivalent width from W3D ~ 0.61 pm to W3D ~ 0.17 pm). 
Moreover, all contribution functions are shifted to slightly higher 
layers because of the presence of saturation effects. The upward 
shift is more pronounced for the 3D contribution function, be- 
cause of the strongly nonlinear fluctuations of the saturation 
factor exp{-Tf }. As a result, S(3D) now becomes smaller than 
S((3D)) and S(1D) in the deepest part of the line-forming re- 
gion. Therefore, the ratio of 3D to ID equivalent width becomes 
smaller than in the weak line limit, where S(3D) > S(1D) over 
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Fig. B.5. Disk-center (p = 1) equivalent width contribution func- 
tions, S(log t c ), of a weak (top) and strong (bottom) Fe n line 
with excitation potential^ = 10 eV, at A 850 nm. The contri- 
bution functions have been computed for a single snapshot of 
the 3D model (solid), the corresponding (3D) average model 
(dashed), and the associated ID LHD model (dotted) used in 
this work. They have been transformed from the monochromatic 
to the Rosseland optical depth scale. In all cases, the saturation 
factor exp{-T( } is properly taken into account (see Eq. lB.U . 



the whole optical depth range. Hence, the total 3D abundance 
correction, A3D-1D, becomes less negative if saturation is taken 
into account. 

If the line strength is increased even more, the contribution 
functions become wider and extend to higher atmospheric lay- 
ers, as shown in the bottom panel of Fig. lBTBI Recalling that the 
equivalent width contribution function is a superposition of the 
line depression contribution functions for the individual wave- 
length positions in the line profile, it seems evident that the 
double peak structure is related to the contributions of the line 
core (left peak) and of the extended line wings (right peak). 
Test calculations confirm this interpretation. Proceeding form 
the top to the bottom of the line formation region, the difference 
S(3D)-S(1D) changes sign from positive to negative to positive 
to negative. Because the contributions of the different layers to 
the abundance correction cancel partially, a straightforward in- 
terpretation of the resulting abundance correction becomes diffi- 
cult. In principle, a detailed analysis of the line depression con- 
tribution functions at individual wavelengths might lead to fur- 
ther insights. Noting, however, that the situation becomes even 
more complicated when considering flux spectra (involving in- 
clined rays), we have some doubts that such an investigation is 
worthwhile. 
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Fig. C.l. Number density of different molecules, normalized to 
the total number density of carbon (sum over all molecules and 
ionization states) as a function of the Rosseland optical depth in 
the ID LHD model. 

Appendix C: Molecule concentrations, line 
opacities, and height of formation 

The equilibrium number density of diatomic molecules with 
constituents A and B, Nab, is given by the Saha-like relation 



Nab = Qab(T) 



Na 



Ne 



U A (T) U B (T) \2nm AB kT 



h 2 



3/2 



. D /kT 



(CI) 



where Na and Nb are the number densities (per unit volume) of 
free neutral atoms (in the ground state) of elements A and B, 
with partition functions Ua and Ub, respectively; the molecule 
is characterized by its mass, yti ab, its par tition function, Qab, and 
its dissociation energy D rcf. lCoxll2000T) . Defining the number 
densities per unit mass as X, = N,/p, where p is the mass density, 
we obtain 



Xab=pQab(T) 



X, 



x R 



U A (T) U b (T) \2nm AB kT 



h 2 



3/2 



Do/kT 



(C2) 



Figure IC.ll shows the number densities of our selection of 
diatomic molecules (normalized to the total number of carbon 
nuclei, Nab I E Nc — Xab/ Z ^c) as a function of Rosseland op- 
tical depth in the ID LHD model used in this work. In the photo- 
sphere (log tr oss < 0), essentially all carbon is locked up in CO. 
The decrease of Xcm, Xnh, Xcn toward lower optical depths is 
a consequence of the density factor p in Eq. dC.2t . The destruc- 
tion of all molecules beyond tr oss a 1 is due to the Boltzmann 
factor exp{Do/kT}; a higher dissociation energy corresponds to 
a steeper drop of the molecule concentration with T. 

The opacity of a line with lower transition level energy x, is 
proportional to XabIQqb expj-^/AT), and thus the line opacity 
per unit mass can be written as 



logK f = logX A +logX B + -log6 + (D -x)0 
+ logp + const. , 



(C3) 



where 6 = 5040/T, and the temperature dependence of the par- 
tition functions Ua and U B has been ignored; the molecular par- 
tition function Qab cancels out. 

If both atoms A and B are majority species (e.g. H, N, O for 
the conditions in our red giant atmosphere), then Xa and Xb are 
constant, and the temperature dependence of the line opacity is 
given by 



dlog K[ 
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(p = const.) or 
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Fig. C.2. Same as Fig. IA.2l but for two weak (artificial) molecu- 
lar lines: a CO line with excitation potential^- = eV (top) and 
a C2 line with x — 4 eV (bottom), both at wavelength A 850 nm. 



^£ =Do _ x+ l^_ (P = const.). (C.5) 
86 A 2 6»lnl0 ' v ' 

In general, the temperature dependence of the molecular line 
opacity is more complicated, because Xa and/or Xb are more or 
less strongly temperature dependent due to ionization and/or for- 
mation of different molecules. In our red giant atmosphere, for 
example, the concentration of carbon atoms is controlled by the 
formation of CO molecules. This leads to a strong increase of Ke 
with temperature {dlog Kg/ 88 < 0) for CH and C2 at tr oss < 1, 
such that these molecules can only form in a narrow region cen- 
tered around log tr oss « (see Fig. IC.2l ). 

Finally, we point out that the molecular lines form in the 
same height range as the lines of neutral atoms and ions. Figure 
IC.2l displavs the contribution functions for the most extreme ex- 
amples. The ground state CO line (top panel) shows the most 
extended formation region, centered around log tr oss « - 1 . The 
contribution function of this line is almost identical to that of the 
ground state Fe 1 line shown in Fig. IA.2l The high-excitation C2 
line (bottom panel) originates from a very narrow formation re- 
gion located in the deep photosphere around log tr oss « 0. The 
contribution functions of the other molecular lines considered in 
this study lie somewhere in between these two extremes; the en- 
tire formation region of the molecular lines is always inside the 
height range covered by our model atmospheres. 



